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An online data processing EIT system 



FIELD 

This Invention relates io a device and methods for electrical Impedance tomography (EIT). and in 
particular to a new construction of hardware and software. 

BRIEF DESCRIPTION 

This patent presents the development of a new Electrical Impedance Tomography (EIT) system for 
online data acquisition and processing of measurements from two-phase flows, whose axial 
velocities can be up to 1 Qms"\ with ±5% precision- To be able to measure two-phase flows at these 
velocities requires data acquisition speeds (and processing) of 1000 frajmes/s per dual-plane in 
order to generate the flow velocity maps. The system has been carefully designed tn a modular 
fashion and can consist of several data acquisition modules and computing modules. The fomner 
incorporates a YoJtage controlled current source with a DC-restore circuit, an over-zero switch 
circuit, an equal-width pulse synthesizer unit to produce synthetic waveforms for electric field 
excitation and demodulation signals (for synchronoation and sampling) and a synchronized digital 
demodulation unit with sixteen parallel measurement channels. The computing module 
incorporates a powerful Digital Signal Processor with an IEEE1394 interface for communication. 
These modules form a data processing pipeline for a series of tasks ranging from flow 
measurement to image reconstruction to determining flow velocities. The DSP system can 
implement three image reconstruction methods, these being the SCG method, the single-step SCG 
Method and the Back-Projection method. An online-updating cross-conelatton algorithm has been 
developed and is to be implemented on the system. This is particularly suited to flow velocity 
calculations, as this technique is less time and memory consuming compared to direct cross- 
correlation techniques. The system architecture, special design techniques, and data processing 
algorithms will be introduced in the paper along with the current progress In the system's 
development 




DETAILED DESCRIPTION 

The invention will now be illustrated by way of example only with reference to the accompanying 
drawings, In which 

Figure 1 illustrates the system architecture; » 
Figure 2 illustrates a diagram showing one possible system configuration of high performance two- 
phase flowmeter using the DSP module as a building block, (a) Single Read-Operation timing, (b) 
Single Write-Operation timing; 

Figure 3 shows the Interfaces of Read-timing and Write-timing where t1* 0-60ns; 12=0-252113, t3-0- 

12ns. The tGS signal in the figure is chip selection signal decoded from address signal EA(17:2]. 

Figure 4 illustrates the architecture of the EIT DASsubsysterns; 

Figure 5 illustrates the structure of Equal-Width Pulse synthesizer (EWPS); 

Figure 6 shows a staircase cosine wave (a) in time domain and (b) In frequency domain; 

Figure 7 shows the programmable phase shift function, (a) phase shift 15°, (b) phase shift 45 , (c) 

phase shift 90°; 

Figure 8 illustrates the structure of Synchronised Digital Demodulation (SDD) 
Figure 9 shows the configuration of the sampling modes, (a) 4 samplings, (b) 8 samplings, (c) 16 
samplings, (d) rolling-sample method; 
Figure 10 shows the reading data timing; 
Figure 11 shows the sampling signal for measurement timing; 
Figured shows an equivalent circuit for the EIT transducer electronics: 
Figure 13 shows a realisation circuit of Over-Cuirent-Zero Switching; 

Figure 14 gives a typical curve of the Independent measurements obtained from a 2D phantom 
filled with water. Data was obtained from a 2D phantom filled with mains tap water; 
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stacked from 1000 EIT images reconstructed using SBP algorithm. The 1000 frame? of data were 
taken in 1.06s, 

The measurement system can consist of up to four DSP modules, two EIT electrode arrays 
connected to two EIT Data Acquisition subsystems (DAS's), and a third DAS containing 
conductivity, pressure, temperature and other auxiliary transducers (as shown In figure 1). 

Control Unit + 

The control unit is a multi-processor system contalntag up to four DSP modules. All DSP modules 
are Identical to reduce the design and ? debugging tasks. The TMS32QC6202B is a fixed-point DSP 
chip with 4 ns instruction cycle time @ 250 MHz clock rate, produced by Texas instruments, and 
has been chosen as the processor for the DSP module (HoJliday, 2003). Any DSP module can 
communicate with a remote PC via an IEEE1394 Interface. This makes the functionality of any 
system very flexible. The IEEE1394 interface is fast and can operate over relatively long-distances. 

The four modules work on a data pipeline processing technique to increase the data-processing 
speed. They 'can- perform different tasks on data in pipeline, e.g, the first module controls data 
acquisition and processes raw data; the second module performs image reconstruction and the 
third module fuses the Image data to obtain flow information. A fourth module Is provided to carry 
out other specified functions. 

One application that th© DSP module Is to be used In Is the development of a dual plane Electrical 
Impedance Tomography (EiTJ measurement system as outlined in Section 1. This system is 
required to measure (i) the distribution of local axial, radial and angular velocities of the disperse 
component, (II) the distributions of the local volume fraction of both components and (iii) the 
volumetric flow rate of the disperse component In two component flows. The system has to be able 
to measure local axial velocities up to 10ms" 1 , local radial velocities up to 2.5ms" 1 and angular 
velocities of up to 45 rada" 1 with an accuracy of It Is expected that the dual plane system will 
be* able to capture 2000 frames/s (1 000 frames per plane). 

Figure 2 shows a block diagram of how the DSP module could be used to Implement such a high 
performance system. DSP1 acts as the system controller arid controls the .'data acquisition process 
of two DAS that are connected to \L The DSP then formate the data and transfer it to DSP2. In 
order to acquire 2000 frames/s, DSP1 has to be able to move 409$ 0000 bytes/s from the DAS' 
using the adjacent electrode data acquisition protocol. The time taken to acquire the 2000 frames te 
difficult to calculate as the DAS' operate at much slower speeds than the DSP and are autonomous 
during the measurement process. The date transferred to DSP2 for image reconstruction (data 
transmitted via McBSPO on DSP1 and received by McBSPI on DSP2) will be only 832 000 bytes/s 
after formatting by DSP1, Besides these tasks DSP1 is also responsible for the IEEE 1334. Tfce 
McBSP can. transfer data at rates up to 12.5Mbyte/s. Therefore if one frame is transferred at a time 
with 100% overhead this would equate to 67|is per frame allowing plenty of time to transfer 2000 
frames In one second. If the data transfer packets were larger then the efficiency of this transfer 
would increase. Where possible the data transfer operations will be carried out by the DSP's direct 
memory access controller (DMA), This is an on-board peripheral that is very similar to the DSP's 
computational processing unit (CPU), but is concerned only with transferring data from one location 
in the memory map to another without the intervention of the CPU. 

DSP2 is charged with Implementing the image reconstruction algorithms on the formatted data 
received from the previous module thus generating image data for each frame acquired by the 
DAS. The time taken for the Image reconstruction for one frame is obtained from the matrix 
multiplication of a [104x104] matrix by a single 104-column matrix where there are 104 
independent measurements and 104 finite. elements In the sensitivity map. This gives 10816 
muitlply-accumulates (MACs). The DSP is capable of 600 million MACs (MMACs) therefore 
10816MACs/50QMMACs results in a time of 21.6ns per frame. If a generous overhead Is added to 
make this up to SOjis/frarne for data accesses then the total time to carry out the Image 
reconstruction for 2000 frames is 100ms. 



The image date generated In DSP2 fs then transferred to DSP3. again via the McBSP. Thfs moduie 
calculates the distribution of the local axial, radial and angular velocities by Implementing a moving 
window of several frames where new Image frame data * arrives end the oldest image date Is 
destroyed. Determining the time required to carry out the velocity calculations depends on the 
actual velocity flow. If the axial flow Is 10ms' 1 this results In only around 3S frames existing in the 
frame window. If there are 3 velocity directions and 104 finite elements tfisn one velocity map will 
take 21.8ns to be calculated. However if we again make this up to 50 us to cover overhead and 
apply this to 2000 frames then we have a total time of 100ms. If the two-phase velocity-flow is 
slower then more frames are used in the velocity calculations. For example if it Is 1ms' 1 then 10 
times mora frames are used. This would take the velocity calculation time up to 1s for 2000 frames. 
This is at the limit of the systems capabilities, however the more than generous overhead should 
allow it to be done. In the system proposed the DSP modules will perform their tasks In parallel on 
a data pipeline. The data acted on in each module is historic In relation to that of the previous 
module as you move from DSP1 to DSP3. 

The DSP External Memory inter-Face (EiWIF) is used to control the DAS's. The related ~EMIF 
signals are: EA[1 7:2]: the low 1? bits of the whole 20-bit address lines, ED[31:DJ: 32 bits data lines, 
/AWE: Asynchronous memory writs-enable, /ARE: Asynchronous memory read-enable, /AOE: 
Asynchronous memory output-enable, TJNP0: Timer 0 or genereUpurpose input, TINP1: Timer 1 or 
general-purpose input. The read-timing and write-tlmlng from the EMIF. as shown in Figure 3, are 
similar to that from the PCI bus of a general PC. This simplifies the modular design and debugging. 
The advantage of the EMIF Is that Its timings are much faster and more flexible due to Ks 
prograrnmability. 

EIT DAS subsystems 

Two identical EIT DAS's are used In this system to collect data from two electrode arrays or a dual- 
plane sensor. A block diagram of She EIT DAS is given in the inside of dash-lfne ©ox (ri Figure 4. 
The DAS'S is composed from four modules: EWPS, SDD, VCCS and a control base. 

An Equal-Width Pulse" Synthesizer unit (EWPS) generates a sine-wave voltage signal, which is 
then converted to two sine-wave currents with a phase difference of 160* by the Voltage Controlled 
Current Source unit <VCCS). The two currents -are used to generate eleettle fields Inside the 
sensor. The Synchronous Digital Demodulation; unit (SDD) measures the electrode voltages end 
converts these analogue signals to digital data.' The digital demodulation design used in this 
system follows that successfully used by earlier Workers (Smith, 1992) that can resolve complex 
Impedances making It suited generally to EIT applications. It then stores all the digital data In a 
FIrst-ln-Firsfr-Out memory (FIFO). The data In the FIFO are than transferred into the DSP modute 
via the DSP Direct Memory Access (DMA) function. TBie SDD end EWPS units are currently 
designed in a modular format which may be integrated into three chips in the future, via ASIC 
design technology, this wilt minimise the systems physical size and also Improve the system 
performance. 

The structure of the EWPS is similar to the widely used Equal-Sample Signal Generator (ESSG). 
An auto-run clock signal triggers the address generator to continuously output addresses to a pre- 
written memory. The memory then sequentially outputs data to a Digital-Analogue Converter 
(DAC). a Jew-pass fitter then smoothes the staircase signal from the DAC to output the required 
sine-wave voltage signal. The difference between the EWPS and the ESSG is that, the former has 
different sample numbers at different frequencies to make the stair-case wave using the same time 
step-length at all frequencies, while the latter has the same number of samples at all frequencies 
so the stair-case wave has different time step-lengths at different frequencies. This difference has 
led to a significant reduction in the complexity of the filter design that the EWPS as it only needs to 
set one comer-frequency for the filter for all signal frequencies while the ESSG needs to set 
different' comer-frequencies for the filter tor different signal frequencies. The specially designed 
EWPS also provides a convenient way to produce a number of synchronized signals that are 
necessary to carry the required control of the SDD. 

16 measurement channels are built into the SDD. Each Channel has one Programmable Gain 
Amplifier (PGA) and one 12-bit Analogue to Digital converter (A/D). The time and phase of the 
measurement sampling of the SDD are funy controlled &y the synchronous signals from the ewps. 
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The A/D conversion and .the data transfer to the FIFO are also controlled by a synchronous signal 
The measurements can be stacked In the FIFO and collected by the DSP module as block data 
once the FIFO fs full. Compared to the conventional sequential and parallel data Collection protocol 
not only does this sub-logical controlled data collection- protocol significantly shorten the 
measurement time and eliminate measurement noise caused by the front-end analogue 
multiplexers, but it also produces a much more accurate sampling time and leaves significant time 
for PSP to perform other processes, , 

A smart AD844 based VCCS with DC restoring techniques Is used In this system, which has a wide 
bandwidth, high output impedance, and excellently balanced large output currents (details will be 
published elsewhere), The two currents are injected into one of 16 pair electrodes via two high- 
performance analogue multiplexers. An over zero switch {OZS) technique has been designed to 
eliminate the transient time caused by the residual charge in the coupling capacitance and 
. electrode-electrolyte interface due to the random switching of the excitation currents. 8y a 
controlling the switching time (the phase difference to the sine wave signal) using the synchronous 
signal with a programmable phase provided from the EWPS. the residual potential after the 
switching is approximately kept is zero. 

To collect one frame of raw data, the time cost by a sequential system and the EIT DASs (based 
on 10 adjacent data acquisition strategy) ran be expressed as Equation 1 and 2, respectively. 

T fr™ = (T^ * + T^c +T fran8 )x N de <N cb -3)/2 (Equation 1) 

(Equation 2) 

where N a | fl : electrode number on one sensing plane, T^ cycle of the exciting sinusoidal current. 
Nsat: cost in terms of cycle numbers for the signal settling, T^qp: acquiring and converting time of 
the A/D converter, T tens : time used to transfer one A/D result into controlling unit and Tf ramfl : time 
used to collect one frame of raw data. 

Assuming T^IOOus {Le. 10kHz exciting frequency). N wl =4 7 Tadc^I^Sus, T^O-OBUS, N ete =16, 
then T tfpm9 of the das will .be 41 ,73, the and 6.42 ms respectively. It can be SBen that the data 
acquisition speed is greatiy slowed down by .ttis.stftling time, especially at low exciting frequencies, 
which is almost linear to the ntimbfer of sfgrM cycles used for overcoming the settling time. The 
application of the novel signal conditioner as described above, the could be set as one and the 
most time of T itwm would be also available for the data processing in the control DSP. which 
significantly improves the data collection speed. ~ 

Signal conditioner 

In previous EIT DAS designs, two digital synthetic (AD700S DOS) chips were used for sine wave 
generation and synchronization (Wang 1B98). One DDS has been adopted to generate a sine 
wave, whilst other generates a sampling signal for demodulation. Although the AD7008 can 
generate a precision sine wave and required sampDng signal, the synchronisation between the sine 
wave and sample signals generated from the two chips are hart to achieve and control. Further It 
is also unable to generate a particular waveform, e.g. a synthetic multi-frequency signal, except a 
sine wave. With the emergence of Impedance spectroscopic examinations of material* a 
simultaneous spectroscopic technique Is becoming desirable (Riu 1995; Shallof 1997). 

Equal-width Pulse Synthesizer (EWPS) is a sine wave voltage generator and synchronous signal 
generator, as shown In Figure 5. This is composed of a programmable look-up table for the sine 
rn»w^^ i? Analogue Converter (DAC), a low pass filter, a 8 bit Digital to Analogue. 

Converter (DAC) and other logic counting circuits. An oscillator and counter an? used to repetitively 
clock the stored 'staircase' wave out of PROM, which is then converted to analogue waveform by 
the DAC and the low pass filter. The amplitude of the sine wave can be controlled by a 0 to 1.25 V 
by an 8 bit DAC. The digital approach can easily generate both the synchronous signal and the 
sample signal. The synchronous signal (square wave) Is used to control erf the siart of the 
measurement. The sample signal fs a trigger for synchronized digital demodulation. 



In the drcurt, a PROM-based staircase function generator has been uitfiseci. The prom i* 
programmed with either a dlgWsed sine wava signal of other particular synthetic waveform 
staircase, In which the wavsfbrm Is stored by an equal width step for all different frequencies. In the 
EWPS operation, the original dock frequency is fixed and unchanging with signal frequency. This 
uses different sample numbers at different frequencies to make the staircase waveform with the 
same step-length of time for ail frequencies. 

It Is known that a cosine staircase In the frequency domain can be presented by Equation 3 {Evans 
1980). 

= Aiy^(JVfe±l)/jrsin/^] (Equations) 

where N is number of the steps to form a kP* frequency and f 0 is the fundamental frequency. 
The time domain and frequency domain are illustrated in Figure e. 



considering the use of the same step-length, the sampling pulse frequency, f 3 , Is unchanged. , 

/, - W or /■ = 4?/- (Equation 4) 



N 

Therefore, 



™-&%[' J % 1 *H£ (Et " Jae ° n5 ' 



Hence, the lowest harmonics from the EWPS Is n fj and the EWPS may use one low pass filter 
with a fixed corner frequency for all waveforms- The step-length used in the EWPS is 0-5 us which 
gives a sampling pulse frequency Of 2.0MHz. 

The synchronous- signal (square wavfe) is used to control of the start and end of the sampling and 
the sampling pulses are triggers-tor talcing a number Of measurements in a specific time sequence 
from an EIT sensor. The^tlmrnff'ahfd width of both signals play a very important rote in the 
achievement of SDD accuracy. Two Sampling modes are created by the EWPS. The Mode 1 Is 4 to 
32 samplings to be selectable per period, and the Mod© 2 is a rolling-sample method that performs 
one sample per period with a fixed phase shift between one and another. 

An over zero switch (OZS) technique has been designed to overcome the transient time caused by 
the residual charge in capacitances in coupling and electrode-electrolyte Interface- In order to 
achieve this function, the multiplexers have to be operated at a precise phase angle. A 
programmable phase shift (between the synchronised signal and sampling pulses) function is 
adopted (Wang 1994). This is achieved by the control of the high 8-bits of the PROM address, 
which selects one of 128 staircase, waves with different pre-calculated phase shifts also stored in 
the PROM. The phase shift can be adjusted from 0° to 90° with a minimum step of 0703° (see 
Figure 7). 

The Synchronised Digital Demodulator (SDD) is a high speed, synchronised parallel digital 
demodulation system, which is composed of 16 sets of Programmable Gain Amplifiers (PGA), a 
12-bR ADC, a Strobed Flrat-ln Rret-out Memory (FIFO) and control logic. Figure 8 shows the 
functional block diagram of the SDD. 

SDD uses the EWPS sampling signals to acquire the measurement data in parallel at precisely 
defined intervals. The Mode 1 can acquire 4-32 samples per period with a precise timing. The 
Mode 2 is particularly designed for sampling at a high frequency due to me speed limitation of the 
conventional A/D converter. The configuration of the sampling modes is shown in Pi9"^- J he 
completion signal of the ADCs with a control clock will then sequentially push the converted data 
into two FIFO memories. All the operations are logically controlled without the involvement of the 
DSP, except the starting signal. 
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Not onJy the parallel data conectton structure can significantly shorten the data collection time but 
also eliminate measurement noisy caused by front-end analogue multiplexers used in the 
sequential data collection structure. The required minimum data acquisition time per sample is 
given as, 

T m « To* +S7^ +7^ (Equation 6) 

where T a la the data acquisition time, dr .* ' • 

Tado is data convert time for AID converter. 
Tr Is the time for reading the data from A/D converter to FIFO. 
Tfctey Is delay time for logic circuit 

Considering ail factors, the minimum date acquisition time is about 2 us, which 1b shown as Figure 
10 

The structure of logical control without the involvement of the DSP Is a lighting point for the signal 
conditioner. After the 'start' command Is set by DSP, SDD will automatically sample the sensor 
outputs at precisely defined intervals within one or more periods and convert the analogue data to 
digital and then, push the digital data into FIFO. The SDD will then return to its initial status (Figure 
11). In order to 'inform' the DSP of the working status, a 'busy* signal and a FIFO memory status 
with the start signal from DSP construct a simple communication interface. This frees up the DSP 
to perform other tasks as well as shortening the acquisition time. 

Conventional, ERT system uses sinusoidal constant current sources to excite the sensing field. 
Since a controlling unit, such as a PC, DSP or microcomputer, directly controls the current 
switching between electrodes, the switching operation may occur randomly during anytime of a 
sinusoidal cycle. When the current value at the tim& Of the Switching Operation is non-zero, 
electrical charge may reside in the capacitances, of the signal passing circuit, as double layer 
capacitance, coupling capacitance and stray capacitance. These may result in residual potentials 
presented between these capacitances and produce an additional contribution to next excitation 
and measurement. These residual charges will require several sinusoidal cycles to disappear, 
which means the whole data acquisition system (DAS) has to wait for several cycles after each, 
switching operation in order to perform next voltage measurement This would oe harmful to a . 
parallel measurement structure. 1 l.f. 

To analyse the effect, an equivalent circuit for one electrode measurement channel is given in 
Figure 12(a). /where C H is the Helmholtz double layer capacitance, R* Hie charge transfer 
resistance, Rb the bulk resistance, C c the coaxial cable capacitance and C© the other stray 
capacitances (Wang 1994). C f the AC coupling capacitance (Bockris 1969; Pollak 1974; Wang 
1994). Taking typlcaJ values of C H of lOpF/cm 2 , Rc of SOQQxm* (Pollak 1974). CfOf IjAFand Rf of 

1 MO. In the majority or cases, Cs and Cc has a value of around a hundred plcoFarads, therefore 
their effect on the solution is minimal and can be ignored without significantly altering the result 

Assuming the current Injection is switched off at to, the voltage at the input of the differential 
amplifier can be expressed as: 

where I a a ^sin(atf 0 + <p) , T x = R 0 C# and T 2 = R f C f . (Equation 8) 
As &R 9 C H » 1 , T t « T 2 , then equation 5a can be simplified to: 



Ki = ~K& c n Y x Ae"*** sin(orf + p) = -(tfC„ ) ^Am** cos(d?/ + <p) (Equation 9) 



A simplified equivalent circuit and fe transient curve are shown in Figure 12(b). Assuming cu«20re 
kHz, y =0. A~10 rnA, and the current is switched off at t=to=0, the transient voltages at the Input of 
the differential amplifier should still be 14.4 mV and 5.89 mV in respect to the settling time At 
=0,1 T n and At =1*! , (Ti -0,005 e). Therefore, (he major contribution to the changed residual voltage 
Is from the double layer capacitance, which results a significant drop of 1 the measurement accuracy. 

The over-zero switching (OZS) discussed here is aimed io eliminate trie charged residua! voltage 
by controlling the switching operation'of the current injection. In an* Ideal case, no charged potential 
resides in either the double layer capacitance or the coupling capacitance after switching off the 

current injection so the settling time Is shorter with the data acquisition speed being increased. 

Since the charged voltage has a u/2 dalay with respect to the current, the switch should operated 
at either zero-voltage between the double-layer capacitance or maximum-purrem through the 
capacitance. The current-maximum is actually applied to optimise the switching time or the relevant 

phase difference between the synchronised signal and the sinusoidal signal, as the charged 

voltage between the double-layer capacitance is hardly measurable. 

Figure 13 is the circuit used in the Improved system to realize the over-eurrent-zero switching 
technique. Two inverse-phased currents C1 & C2 are connected to the common Inputs of two 
multiplexers U1 & U2 respectively. 16 electrodes mounted in one sensing plane are connected to 
the 16 output channels of each multiplexer. Channel selection is together controlled by DSP 
command and a synchronize signal, SYN.SIN, of current C1 . SYN_SIN Is a square-wave signal, 
which has the same over-zero point with two current signals C1 & C2. When DSP send out a 
switching command, an $-bits data (each multipbcer is controlled by 4-blts data), EDtf):7], is written 
into the first flip-flop, U3. This data is later written into the second flip-fiop, U4, by. the next rising 
edge of SYNJS1N. and finally make the two muitlplxera switch to other channels. A D-type flip-flow, 
USA. is used In this circuit to prevent writing data into U4 before the data has been set-up at the 
output of U3, 

Data processing algorithms 

Three kinds of image reconstruction algorithms have been coded for the DSP module.. They are: 
SCG method* (Wang, 2002) f Sensitiviry-map-basad Back Projection (SBP) method (Kbtre, 1904) 
■* :y ancfca simpjifi&d. SCG method - Single-step SCG (SSCG) method. The SCG method .b-.a-fclci^gj ./ : .. r . * 
liej&tive method, which has the highest-performance but Is the most complicated' end tnHS- " - 
consuming of the three. The SSCG method is the SCG method iterated only once from a 
homogenous conductivity distribution. The SBP method is the simplest and fastest; It seems that 
the SSCG has a better performance than the SBP method. 

The time taken to reconstruct one frame Image will be different for each of these algorithms. As an 
example, the time taken using the SBP reconstruction for one frame is Obtained from the matrix 
multiplication of {104x104] matrix by a single 104-column vector where there are 104 independent 
measurements and 104 finite elements in the sensltMty map. This gives 10816 muttiply- 
accumuiates (MAC). Based on the DSP's capacity of 500 million MACs per second, it needs 21.6 
us to complete the matrix multiplication. Adding 100% overhead time, a generous image 
reconstruction rate of 50us/frame has been calculated. Then the total time to carry out the Image 
reconstruction for 1000 pair of Images (2000 frames) Would be 100 ms. 



Assuming that the two EIT electrode arrays are installed parallel on the flow^pipe, the axial-flow 
velocity distribution can be estimated by a direct cross-correlation method, as given in Equation 1G: 

HaW- E/i<"-»>/2<<*>' -.W-l> (Equation 10) 

where N Is the sample length, and frfm) & f 2 (m) are the m* up-flow and down-flow images 
respectively. Equation 10 can be simply implemented online by updating the R12M with the new k- 
th images, as described In the following: 

^W»^W)+/i»-«)/aC«- n»o,i,~-,(N-t> (Equation 11) 
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when* superscript (k) indicates that the corresponding value Is of K-th sample. This Implementation 
can greatly save calculation lime and reduce memory size. The Compulation speed fe subject to the 
number of pairs of Images that ere estimated as a maximum sampling length- For example fake a 
velocity ranging from 1 m/s to 10 nVs giving a turndown ratio of 10:1 and V max *1Qm/s, V mtn =1m/s. 
Assuming the distance, D, between the two sensing planes In a dual-plane sensor Is 0,1 m and the 
data acquisition speed, V**,. is 1000 frames /sea the sampling length, N, Is 2a frames far and 
200 frames for V^! The minimum number of multiplications required to find k pixel flow velocity in 
the velocity range is . 



MAC, 




* * 1 80 (MAC) , (Equation 12) 



which give 37.44^s for generating en axial velocity map with 104 pfcels and SO jis if some 
overheads are expected. To secure the data accuracy, a repeated computation is normally applied 

to the sampling window. For a total number of 2000 frames data, it needs 100 ms for the 
completion of the axial velocity map. According to other workers (Beck, 1986; Wah, 1998), the 

minimum acquisition time, 5 T can be taken as twice the product of the minimum transit time of the 
fluid, t, and the fractional velocity discrimination, tc. By rearranging gives the following 

k = 6/2t (Equation 13) 

If the distance between two sensing planes is 0.1m then x ° 0,1m/10rns M = 0.01e, If 1000 frames/s 
was achieved then 5 - 0.001s. Using Equation 6 gives a discriminatory precision of 0.05 or 5%. 

Preliminary result 

The fast data collection and implementation system has been designed end built at the university of 
Leeds. A number of trials have been Successfully performed. Without making efforts to Optimise the 
control code for data collection, transferring and processing, a data collection rate of 1000 frames 
per 1.06s was recorded. A 2D phantom fitted with 16 Stainless Steel Electrodes : and filled with 
• -mains tap water was used in one of the trials. In the experiment, all this PQXsffifcB SDBtiSre 
equally set to a fixed gain value In order to Simplify the test The Current excitation frequency was 
80 kHz. Measurements were taken from all electrodes during each current exdtettbnl^rojection, 
which gave 256 measurements. Four samples per period were taken. Therefore, each frame data 
has a total of 1024 measurements. A simplified digital matched-filter was applied to demodulate the 
data. A set of Independent measurements were then extracted and averaged from their mirror data 
based on the reciprocal theorem. The new EIT System has not been optimised at the time of 
carrying the initial test. It is expected that a rate of 1000 pairs (2000 frames) per second could be 
achieved when two DAS are connected to the system. 

The top and middle curves In Figure 14 present two measurements obtained from a vessel filled 
with mains tap water. The bottom curve denotes the relative change between the two 

measurements to reflect the repeatability of these measurements. Figure 15 shows an Image of a 
high conductive (100mS) brine droplet (approximately 5ml) injected into a mains tap water 

background. The Image was stacked from 1000 2D images and the droplet was visualised from the 
axial cutting section of the stacked image. As this Is preliminary work, validation of these 
measurements will be made In the near future and reported on. 
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Figure 1: System Block Diagram 

Figure |: A diagram showing one possible system configuration, of «]high 
performance two-phase flowmeter using the DSP module as a building 
block, (a) Single Read-Operation Timing, (b) Single Wrlte-Operatlon 
Timing 

Figure 3: Interface Read-Timing And Write-Timing Where t1= 0-60ns; t2=0-252ns, 
t3=0-12ns. The /CS signal in Fig .2 is chip selection signal decoded from 
address signal EA[17:2], 

Figure 4: Architecture of EIT DAS 

Figure 5: The structure of Equal-Width Pulse Synthesizer (EWPS); 

Figure 6: A staircase cosine wave (a) In time domain and (b) in frequency domain 

Figure 7: Programmable phase shift function, (a) phase shift 15°, (b) phase shift 
45°, (c) phase shift 90°; 

Figure 8: The structure of Synchronised Digital Demodulation <SDD) 

Figure 9: The configuration of the sampling modes, (a) 4 samplings, (b) 8 
samplings, (c) 16 samplings, (d) rolling-sample method 

Figure 10: Reading data timing 

Figure 1 1 : Sampling signal for measurement timing 

FIgure12: An equivalent circuit for the EIT transdliqsr electronics 

Figure 13: A typical curve of the Independent measurements obtained .from a 2D phantom filled 
with water. Data was obtained from a 2D phantom filled with mains tap water. 

Figure 14: An axial cutting image of a high conductive droplet dropped Into water, 
which Is stacked from 1000 EIT Images reconstructed using SBP 
algorithm. The 1000 frames of data were taken in 1 -06s 
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Figure 1: System Block Diagram 



Inter-processor communication 




Figure 2: A diagram showing one possible system configuration of high performance two-phase flowmeter using 

(he DSP module as a building block. 
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(b) Single Write-Operatlon Timing 



Figure 3: Interface Re^Tlmlng Antf Wrlte-Tlmlng Whore tf» O-Wne; tZ=0-262ns, t3=0-12ns. The /CS signal In Flg.2 
Is chip selection signal decocted from address signal EA[17^2]. 
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Figure 5: Tha structure of Equal-Width Pulse Synthesizer (EWFS) 
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(a) (a) 
Figure 6: A staircase cosine wave (a) In time domain and (b) In frequency domain 
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Figure 7: Programmable phase shift function (a) phase shift 1 5° (b) phase shift 49° (c) phase shift 90" 
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Figure B: The structure of Synchronised Digital Demodulation (SDD) 
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Figure 9: The configuration of the sampling modes (a) 4 samplings; (b) 0 Samplings; (c) 16 samplings. 

(d) roJIlng-sampIo method. 
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Figure 11 : Sampling signal for measurement timltig 





Flgurel 2: An equivalent circuit for th* ETT transducer electronic 
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Figure 13: Realisation Circuit of Over-C*lrrent-Zero Switching 
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Figure 14: A typical cutVo of the Independent measurements obtained from a 2D phantom filled with water. Data 
was obtained from a 2D phantom flllod with mains tap wgfer- 
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Figure 1 5: An aadai cutting Image of a high conductive dropta dropped Into water, which is stacked from 1000 ETT 
images reconstructed using SBP algorithm. The 1000 frames of data were taken In 1*05 ft 
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